Magnetization suppression of Type-II Superconductors by external alternating 

magnetic field 
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The effect of suppression of static magnetization of an anisotropic hard superconductor by al- 
ternating magnetic field is analyzed theoretically. The magnetic moment suppression dynamics is 
described with respect to the magnetization loop of the superconductor. It is found that in some 
cases the magnetic moment varies nonmonotonically with the growth in amplitude h of the alter- 
nating field. Effect of transition, induced by h(t), of superconductor form paramagnetic into the 
diamagnetic state is considered. The amplitude of alternating magnetic field hc{S,'&) for which 
the complete suppression of the magnetization occurs is calculated as the function of anisotropy 
parameter 5 and it orientation angle i9 with respect to the crystallographic axes of the sample. 

PACS numbers: 74.25.Ha, 74.25.Sv, 74.25.Qt 



I. INTRODUCTION 

For many years the physics of the vortex matter in 
the superconductors attracts the attention of researchers. 
This is because collective properties of the vortices reveal 
variety of reach and very interesting phenomena which in- 
volve different phases and phase transitions RefiSiSi^, 
magnetic instabilities Ref.^-^'^'^ mesoscopic and fluctua- 
tion effects and others. In this paper one example of the 
instability effects, namely the problem of magnetic mo- 
ment suppression by external alternating magnetic filed, 
is considered in some details. Usually the static and 
quasi-static electromagnetic properties of hard supercon- 
ductors is treated in terms of critical state model at first 
suggested by Bean-'^". This model describes of magnetic 
induction B distribution inside the superconductors. In 
the simplest possible case, when superconductor repre- 
sents rectangular slab with x-axis directed perpendicular 
to it planes and z-axes along the magnetic field direction, 
the Bean equation has following form dB/dx — zLAttJc/c, 
here Jc is critical current density. Critical state equation 
can be easily understood by analyzing the forces balance 
which acts on a single vortex. In accordance to the model, 
the Lorentz force, which acts on the vortex, is compen- 
sated by the pinning forces between the vortex and vari- 
ous defects in the crystalline lattice. Later critical state 
model was generalized for more complicated situations in 
several directions. 

First generalization was performed by the authors of 
Refjii in order to include quasi-stationary processes. In 
the most of the cases, for simple geometries, direction of 
the critical current vector is uniquely determined. But 
situation is some what more complicated when external 
magnetic field has several spatial components or when its 
varies in time. For these time dependent cases usually 
critical state model is written as curlB = (47rJc/c)E/i?, 
where E is electric field. This model assumes that vary- 
ing in time magnetic flux produces electric field, super- 
conductor becomes in the resistive state and the direction 
of electric current coincide with the direction of E similar 
to normal metal. It is also important to answer the ques- 



tion about the direction of critical current density in that 
regions of superconductor where there is no electric field. 
At this point rises one of the most important electrody- 
namic properties of hard superconductors. The magnetic 
state of the superconducting sample, at some specific mo- 
ment of the time, is described not only by critical model 
equations with some boundary conditions but also de- 
pends on prehistory - how this state was prepared. In 
order to understand this statement lets consider example 
of static magnetization. Assume that we have supercon- 
ducting slab placed in the external magnetic field with 
the quasi-stationary varied amplitude. Because magnetic 
field changes in time very slowly we can assume that in 
each instance superconducting sample is in the critical 
state with the critical current density Jc(-B) which is 
determined by the value of magnetic induction at this 
moment of the time. Direction of this current coincides 
with the direction of electric field E which emerges from 
slowly varying external magnetic field H. This picture is 
valid while amplitude of H changes. As soon as magnetic 
field stop changes the electric field disappears but critical 
current density persists its direction which was set by the 
electric field at last instance of its existence. Exactly this 
dissipationless current determines final magnetization of 
the sample which is usually called static. 

It is also important to point out that critical state 
model equations are significantly nonlinear and this non- 
linearity is peculiar only for superconductors and has no 
equivalents in other nonlinear medias. This specific non- 
linearity leads to several very interesting cffectsAi one of 
which is static magnetization suppression. The essence 
of this effect is the following. Let a plane superconduct- 
ing sample cooled in zero magnetic field be placed into 
an external magnetic field H > Hd (here Hd is the 
lower critical magnetic field) which is parallel to the su- 
perconductor surface. Pinning leads to the emergence 
of a nonuniform distribution of the magnetic induction 
and, accordingly, a static magnetization, in the sample. 
If an alternating magnetic field h.{t) — hcos(tjt) is ap- 
plied to the magnetized sample in a direction parallel to 
the sample surface and perpendicular to a constant field 
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H, then magnetization M of the sample decreases. Di- 
rect measurementsiSi show that everywhere in the sample 
where the alternating magnetic field penetrates, the flow 
of nondissipative currents becomes impossible. The cur- 
rent that previously screened field H and contributed to 
the magnetization of superconductor disappear. If the 
amplitude is sufficiently large h ^ Hp = 2ttJc{H) /cd 
(the alternating field is practically penetrate through the 
entire sample) magnetization is suppressed. 

Second generalization of critical state model is related 
to the fact of high anisotropy of the high-Tc supercon- 
ductors which will sufficiently important in this work. In 
the anisotropic model critical current density becomes a 
tensor Jdj and straightforward generalization of model 
equations gives (curli?)^ = {4Tr/c)JcijEj/E. Diagonal 
components of the Jdj tensor are significantly different, 
for example for YBaCuO the critical current density in 
the ab plane is much grater the that in c direction. 

In this work the theoretical analysis of static magne- 
tization suppression in the anisotropic case under influ- 
ence of alternating magnetic field h(i) — hcos(ijji) is per- 
formed. It is known that magnetization curve M{H) of 
HTSC reveals hysteresis and it turns out that character 
of magnetization suppression is very sensitive to the po- 
sition of current magnetization in the magnetization loop 
(this reflection of the sensitivity to the magnetic prehis- 
tory described above). It is found that in total there are 
nine regions in the magnetization loop where suppression 
scenario is qualitatively different and dynamics of mag- 
netic moment suppression with growth of the amplitude 
h is described in each of these regions. 



II. MAIN EQUATIONS AND GEOMETRY OF 
THE PROBLEM 



Consider an infinite plane-parallel superconducting 
plate of thickness d placed in external constant magnetic 
field H and alternating magnetic field h.{t) = hcos(a;i) 
which are mutually perpendicular and parallel to the 
plate surface. It is assumed that all fields and currents 
depend on only one spatial coordinate x directed along 
the normal to the plate. The origin a; = is located in 
the origin of the sample. In this geometry the equations 
of generalized critical state model for magnetic induction 
B, written in the components, take form 



dx 



in 



Jcy{By,B^)cOs(if{x)), 



dBy Att 

= —Jcz{By,B^)sm((p{x)), 



(la) 



(lb) 



here f{x) - angle between the electric field vector E and 
y axis. The spatially averaged component of the magne- 
tization M along the direction of the external magnetic 
field H, which we will denote as Mh, is given by following 



formula 



Mh = — 



Aw 



1 fd/2 

- / {B,{x)coa^ + By{x)sim})dx- H 

" J-d/2 



(2) 

Eqs. and ^ should be accompanied by the Maxwell 
equations for electric field 



dx c dt 
and boundary conditions 

d' 



dEy 

dx 



c dt 



(3) 



B, ± 



Hcos-d. B„[± 



Hsin'd. (4) 



It is clear that exact integration of Eqs. is very compli- 
cated problem due to specific nonlinearity of equations. 
For the sake of simplifying of the calculations, the depen- 
dence of the components Jc(B(x)) of the critical current 
density on the x coordinate, caused by the nonuniformity 
of the magnetic induction distribution, will be neglected 
and assumed that Jc (B(a;)) = Jc(H). In addition, the x 
- coordinate dependent angle ^{x) between the electric 
field and the y axis will be replaced by the angle TT/2 + 'd 
between the external alternating magnetic field and the 
z axis (as in the anisotropic situation). It turns out that 
these simplifications have no qualitative influence on the 
results but make it possible to perform analytical calcu- 
lations completely. 

For further convenience lets introduce following dimen- 
sionless variables 



B 



B 

Hp ' 



Hr, 



Hr, ' 



H„ ' 



M - 



Mh 

Hr, 



indJr 



(5) 



here capital letters refer to the dc-magnetic field H and 
small letters to the ac-magnetic field h{t). 



III. CALCULATION OF MAGNETIC MOMENT 
M SUPPRESSION DYNAMICS 

The main idea of calculations is as follows, first of all 
we have to find static magnetization when the alternating 
magnetic field is absent. This part of the problem is well 
developed both theoretically and experimentally and we 
can refer to Reftii where magnetization loop was studied 
in details. For each point of the magnetization loop there 
is some specific distribution of magnetic induction which 
determines magnetization. As the second step we have 
to solve dynamical problem and describe the penetration 
of ac-magnetic filed. In this way it will be possible to 
find the penetration depth and critical amplitude of 
ac-field /ic(^) a-t which total suppression of magnetization 
occurs. As soon as problem is two-dimensional it means 
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that from two penetration-depths of each spatial compo- 
nent of alternating magnetic field we have to choose the 
highest. 

One can easy shows that at all places where the alter- 
nating field h{t) penetrates, the magnetic induction can 
be presented as a sum of two terms. One of these terms is 
a constant homogeneous quantity coinciding with vector 
Ti. The second term, which corresponds to the nonhomo- 
geneous magnetic induction distribution of ac-magnetic 
field by^z, may be described by the following equations 
(in analogy with Eqs.lQl) 



— sm — = 6 cos 'd. 
04 ot, 



(6) 



These equations hold in the superconductor region where 
both induction components bz and by are present. In the 
region there component bz vanishes, which correspond to 
ip{x) = 7r/2 or equivalently = 0, and there is only in- 
duction component by, the distribution of this component 
is described by the equation 



dby^ 

Solution of Eqs.JSJ and Eg.l!?!) is easy to find 



(7) 



bziO ^hsin^ + sini}{^~l), 6 < ^ < 1, (8a) 



r /icosi?-f (5cosi9(?-l) <C<1, 

'^^^^^"l byi^z) + - ^z) ^^^> 

here £,z and ^y are penetration depths of each compo- 
nent of ac-magnetic field. Cusp in the distribution of by 
component is the result of the fact that coordinate depen- 
dence of angle (p{x) was neglected. From the condition 
bz{S.z) — one finds the penetration depth for z compo- 
nent of ac-magnetic field ^z — ^ — h and value of by at that 
point by{S,z) = (1 ~ S)hcos'&. Similarly one can find pen- 
etration depth for y-component S^y = 1 — /i— {-^-j^) ft. cost? 
and finally from the condition 6j^(^ = 0) = critical am- 
plitude /ic(i^) of the alternating magnetic field at which 
total suppression of magnetization occurs 



(l-d) cost? 



(9) 



Function hd^) has universal character; it is defined via 
the anisotropy parameter 5 of the theory which makes 
this formula to be interesting from the point of view of 
experiment. Say by measuring the hc{i)) a,t •& — Q one 
can directly find anisotropy parameter because 5 = he or 
5 = chc/2'KdJcy in dimension variables. 

Let now consider the dynamics of suppression of the 
static magnetic moment M. of the sample by an orthogo- 
nal alternating magnetic field. As has been noted above, 
the suppression effect essentially depends on the mag- 
netic prehistory of the sample, i.e. on the position of the 



starting point on the magnetization loop. Below will be 
considered the most simple but nevertheless practically 
interesting case of the magnetic prehistory, when the ex- 
ternal magnetic field monotonically increased up to cer- 
tain maximal value Tim, such that 1 <C Ti-m ^ '^c2, and 
then decreased back to zero. We will distinguish between 
two cases namely direct magnetization, when TL monoton- 
ically increased Ti G [0, Tim] and reverse magnetization 
when magnetic field Ti. gradually decreased from its max- 
imum value to zero TL € \Hm , 0] . It turns out that there 
are three distinct regions for the direct magnetization and 
six for the reverse where suppression occurs qualitatively 
and quantitatively different all these regions will be dis- 
cussed. For each of these regions the dependence of the 
magnetization on the amplitude of alternating magnetic 
field M{h) will be found. 



A. Direct Magnetization 

For the direct magnetization there are three distinct 
regions H < hc{d), h^'d) <H<landl<n< Hm- 
The calculations of A4{h) are straightforward but cum- 
bersome so that the only final result will presented and 
discussed. 

A.l Range: < H < hd^) 



0<h< hc{^)~U 



2£_ 

87T 



8irft,c(i?) 



(10a 



hc{i})~n< h < hd^){i-n), 



8tt 



dH? sin 5 sin 

87r/ie(i?) 



SHsmi) _ 2i. 



M{h) = -gL(cos2 ^ + Ssin' (l - 4^)^ 



(10b) 



(10c) 



A. 2 Range: ft.c('?) < < 1 - if amplitude of the ac- 
magnetic filed is at the range < h < /ic(t?)(l — H) 
then magnetization M(h) is described by the formula 
Ea. (fTOb|l . For a higher amplitudes hci-d){l - H) < h < 
hc{i}) the M{h) is given by Ea.lfTnc|l. 

A. 3 Range: 1 < Ti. < Tim ^ for any amplitude of the 
alternating magnetic field at the interval < /i < hc{'d) 
the suppression of M.{h) occurs in accordance with the 
formula Eq. (|10c|) . 

Generally speaking the dynamical suppression of Ai 
for the direct magnetization occurs by the same sce- 
nario for all three regions in the magnetization loop 
TL e [0,/ic(t?)], TL e [ftc(i^),l] and TL e [1,H,„]. Because 
of that we will discuss only the case A.l and for two oth- 
ers situation is essentially the same. Inside of the region 
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Ti G [0, /ic(^)] we can distinguish three steps in the evo- 
lution of the magnetic induction distribution, and subse- 
quently magnetization in the sample. Each of these 
steps occurs at some specific range of amplitudes of alter- 
nating magnetic field which is described by inequalities 
in the formulas Eas. Hl()a|l - Hl()c|l . Static magnetization is 
negative which means that sample is in the diamagnetic 
state. With the gradual increase of h the static mag- 
netization is smoothly vanishes. This happens because 
ac-field penetrates dipper inside the sample, magnetic in- 
duction in that region becomes homogeneous and doesn't 
contribute to the magnetization. Full suppression occurs 
at the critical filed Eq.©. 



B. Reverse Magnetization 

Suppression of A4 for the reverse magnetization is qual- 
itatively different from that of direct magnetization. For 
all magnetic filed ranges in the magnetization loop (only 
with exception for the last one) there is characteristic 
cusp in the distribution of the magnetic induction. This 
cusp is the result of redistribution the magnetic induc- 
tion after lowering amplitude of H below the TLm and 
there is peculiar feature of critical state model. Precisely 
this cusp is responsible for some new effects. This hap- 
pens because now static magnetization has two terms 
one of which is negative but another is positive. De- 
pending on the value of 7i and h one of this terms wins 
such that magnetization may be either positive (para- 
magnetic) or negative (diamagnetic) and paramagnetic- 
diamagnetic transition is possible. 

B.l Range: Hm - S < H < Hm 

0<h<he{^)iH,n-n), 

M{h) = -^{cos^i) + Ssm'd) (l-^)-f 

(11a) 



-H 



4:TT 



WttS 



M(h) 



n) <h<s-^h,{i)){n,n 

^ 2/ic(i5) 



47r 



n), 



16n5 



Stt 



1 



(lib) 

In the amplitude interval S ^hc{i!)){Hm — 'H) < h < hc{'d) 
the magnetization is described by the formula Eq. (|10cll . 

In the case of B.l we again can distinguish three steps 
in the magnetization dynamics. Despite the existence of 
the positive term in the static magnetization the sample 
still diamagnetic. As amplitude of alternating field in- 
creases sample becomes even more diamagnetic because 
at the beginning the suppression affects only positive 
component of the magnetization. After positive compo- 
nent is suppressed completely, magnetization reaches its 
minimum negative value. Further suppression occurs by 
the same scenario as it was for A.l Ea. (|10c|) . 



B.2 Range: - 2(5/(1 + 5) < H < Hm - In 
the interval of the amplitudes < /i < hc{'d){'Hm — 'H) 
magnetization is defined by the formula Ea. Hlla(l and in 
the interval K[d){n„^ - H) < h < h^{'d){2 - S-\H^ - 
H)) by the formula Ea. l)llb|l . For the final interval we 
have 



hei^){2 - s-\n,n -n))<h< hci-d), 



M{h) = -^(cos2?9-,5sin2,9) f 1 - 7^ 



(12) 



Here dynamical picture of magnetic moment suppression 
is the same as in B.l. 

B.3 Range: Hm - '25 < H < Hm - 26/(1 + 6). In the 
amplitude interval 0<h< K{'d){2 - 6-^{Hm - H)) the 
M.{h) is described by the formula Ea. llla|l and for the 
next interval we have 

K{d){2 - 6-\Hm -H)<h< hci^){Hm - H), 



M{h) 



1 



47T 

h 



cos^ ^ 



167r 



8 IT 



1 



■2h'i(-d) 



(13) 

H) <h< hc{d) 



And finally for the interval hc{'d){Hn 
is described by Ea. ((T^ . 

B.4 Range: - 1 < 7^ < - 2(5. For the case < 
h < hc{^){Hm - H) for M{h) we have Eq.lO and for 
the interval hc{'d){Hm — H) < h < hd'd) magnetization 
is described by the formula Ea. ((T^ . 

As soon as physical processes are similar for the ranges 
B.3 and B.4 we will discuss them together. Because of 
relatively big contribution from the positive component 
of the magnetic moment, sample at the beginning is para- 
magnetic. Increase in the amplitude of the alternating 
magnetic field causes decrease in the magnetic moment 
because at the beginning suppression affects only positive 
component of the magnetization. At some characteristic 
amplitude of ac-filed, which is smaller then hc{'d), magne- 
tization becomes zero. This situation correspond to the 
case when contributions from positive and negative com- 
ponents to the total magnetization are equal. For further 
increase of the amplitude h sample becomes diamagnetic, 
magnetic moment reaches its minimal value then growth 
back and suppress completely at the field hc{'&)- 

B.5 Range: - 2 < H < - 1. For the am- 
plitudes < h < hc['d){2 — [Hm — H)) magnetiza- 
tion is described by Ea. H13|) and for the final interval 
hc{'&){2 — {Hm —H)) < h < hci'd) we have new formula 
for the the magnetic moment M{h) 



h,{d){2-{Hm-H)) <h<h,{d), 



M{h) = 



_ cos'' i?+i5 sin'' 1? 



(14) 



B.6 Range: < H < Hm — 2. For all possible ampli- 
tudes of alternating magnetic field < h < hc{d) mag- 
netization A4{h) s described by the Eg. 1)1411 . 

Behavior of Ai{h) in the ranges B.5 and B.6 is com- 
pletely different form that considered above. For these 
ranges there is no cusp in the distribution of magnetic 
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induction components in the sample and magnetization 
is always positive for any interval of h so that sample is 
totally paramagnetic. Growth of the amplitude of alter- 
nating field smoothly suppress magnetic moment which 
nevertheless remains positive. 

IV. CONCLUSIONS 

Theoretical analysis and recent experimentsSiSii^ 

shows that the alternating magnetic field exerts sig- 
nificant influence on the static magnetic properties of 
anisotropic disordered high- Tc superconductor. Switch- 
ing on a sufficiently strong field h(i) orthogonal to the 
static magnetizing field results in the complete suppres- 
sion of the magnetic moment of the sample. The rea- 
son for this suppression is that, at all places where the 
alternating field penetrates, levelling of the distribution 
profile of the static magnetic induction is observed. In 
the other words, in the same spatial region of the sam- 
ple, the constant and alternating screening currents can- 
not coexist. In the conditions when the alternating field 
penetrates into the entire volume of the sample, com- 
plete suppression of the static magnetization takes place. 



The nature of the magnetization suppression, which con- 
sists in a local effect of the mutual influence of differ- 
ent components of the critical current density vector, is 
manifested in the anisotropic situation in a rather pe- 
culiar way. Different components of the magnetic field 
penetrate at different depths, since the are screened by 
critical current densities components of quite different 
magnitudes. This is the reason why the magnetization 
suppression is primarily caused by the alternating filed 
component deeply penetrating into the sample. As the 
result, the anisotropy induces a quite interesting effect: 
to suppress large magnetic moment a small amplitude 
of the alternating signal is sufficient. In the paper, the 
dynamics of the magnetization suppression with increase 
in the amplitude of the alternating filed h is studied in 
details and results are in agreement with that obtained 
in the experiments Ref^^^ii^. It was shown that in some 
cases the dependence of the moment on h is nonmono- 
tonic and, in addition, during the suppression transition 
of the sample from the paramagnetic state into diamag- 
netic state sometimes occurs. All results can be inter- 
preted within the framework of a critical state model 
generalized to the anisotropic case. 
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